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ABSTRACT 

Model stellar atmospheres are fundamental tools for understanding stellar observations from interferometry, microlensing, eclipsing 
binaries and planetary transits. However, the calculations also include assumptions, such as the geometry of the model. We use 
intensity profiles computed for both plane-parallel and spherically symmetric model atmospheres to determine fitting coefficients 
in the BVRIHK, CoRot and Kepler wavebands for limb darkening using several different fitting laws, for gravity-darkening and for 
interferometric angular diameter corrections. Comparing predicted variables for each geometry, we find that the spherically symmetric 
model geometry leads to different predictions for surface gravities log g < 3. In particular, the most commonly used limb-darkening 
laws produce poor fits to the intensity profiles of spherically symmetric model atmospheres, which indicates the need for more 
sophisticated laws. Angular diameter corrections for spherically symmetric models range from 0.67 to 1, compared to the much 
smaller range from 0.95 to 1 for plane-parallel models. 
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r$ 1. Introduction 

Stellar limb darkening is an important tool for interpreting in- 
terferometric, microlensing and eclipsing binary observations of 
red giant and supergiant stars. It also provides critical infor- 
mation about the te mperature structure of a stellar atmosphere 
, (ISchwarzschildlll906l) as well as a measure of the radial exten- 
sion of an atmosphere dNeilson & Lesterll2012l) . 

Interferometric observations measure the angular diameter 
of a star as well as the intensity variation across the stellar sur- 
face. Some of the first interferometric observations measured 
only uniform-disk angular diameters, that is the angular di- 
amet er for a star assumed to have a constant surface b right- 
i ness (lHanburv Brown et al.l[l974 . IWittkowski et all (12004 pre- 
sented /iT-band interferometric observations of the M3 giant \fi 
Phoenicis with measurements of the first and second lobes of the 
visibility curve, which constrain limb darkening. Unfortunately, 
these observations were not precise enough to distinguish be- 
tween different model stellar atmospheres. Advances in interfer- 
ometric observations ha ve allowed for obse rvations of convec- 
tive cells in Betelgeuse (Haubois et al. 2009) and m easurements 
of gravity darkening in Alt air (Ivan Belle et al] 2001 1 ) . In terms of 
model stellar atmospheres, Aufdenberg et al. ( 20051) constrained 
three-dimensional models using observations of Procyon. 

Microlensing observations, like interferometry, also probe 
stellar limb darkening, but unlike interferometry, which tar- 
gets specific nearby stars, microlensing observations are ran- 
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dom. lAn etall (120021) and lFields et al.l d2003l) constrained non- 
linear limb-darkening relations from microlensing observations 
of a K3 giant and compared them to model stellar atmo- 
spheres. They found significant disagreement between the ob- 
served and predicted limb darkening relation. More recently, 
however, microlensing observations have on ly constrained linear 
limb-darkening relations for red giant stars (Fouque et al. 2010; 



limb-darkening r 
IZub et alJl201lh . 



* Tables 2 -17 are only available in electronic form at 
the CDS via anonymous ftp to cdsarc-u-strasbg.fr or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+ A/ 



Eclipsing binaries and planetary transits provide yet another 
avenue for measuring stellar limb darkening. In terms of red gi- 
ant stars, there are a number of known eclipsing binary systems, 
specifically the £ Aurigae systems that have a K4-5 red giant 
primar y and a main-sequence B-type companion. lEaton et all 
(2008]) fit the orbits for several of these systems assuming a 
simple linear limb-darkening law. There is also the potential of 
observing planets transiting red giant stars, which would pro- 
vide powerful constraints of theories of planetary evolution. 
Currently, extrasola r planets have bee n observed orbiting dwarf 
and subgiant stars dHowell et al but n ot giant star s; future 

missions such as PLATO may remedy this (iCatala et al.ll2010l) . 

These three types of observations are ideal tools for prob- 
ing stellar atmospheres and constraining the physics employed in 
numerical models. Likewise, predictions from model stellar at- 
mospheres h e lp constra i n these types of observations. Recently , 
ISing) d2010h . iHowarthl d2011al) and IClaret & Bloemenl (1201 lh 
presented limb-darkening laws fit to plane-parallel model stellar 
atmosphere i ntensity profiles. Even more recently, IClaret etail 
(I2012L 120131) fit limb-darkening laws to spherically-symmetric 
PHOENIX model stellar atmospheres of cool brown dwarf stars. 
In this work, we study how the assumed geometry of the model 
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stellar atmosphere, plane parallel versus spherically symmetric, 
affects predictions of stellar limb darkening, gravity darkening 
and interferometric angular diameter corrections. We examine 
model atmospheres spanning the effective temperature and grav- 
ity range consistent with yellow and red giant and supergiant 
stars. Tables of limb-darkening and gravity-darkening coeffi- 
cients, as well as new angular diameter corrections are presented 
as more physically based tools for understanding these bright 
stars. 

In Sect. 2 we describe the stellar atmosphere code used in 
this work, as well as the model atmosphere grids computed 
for both plane-parallel and spherically symmetric geometries. 
In Sect. 3 limb-darkening coefficients are presented for several 
commonly used limb-darkening relations. We compute gravity- 
darkening coefficients in Sect. 4 and angular diameter correc- 
tions in Sect. 5. Computations in these three sections provide 
insight into how intensity profiles depend on the assumed model 
geometry that can be directly compared to observations. 

2. Model Stellar Atmospheres 

Model stellar atmospheres form a key foundation of our un- 
derstanding of stars, arguably a great success of computational 
astrophysics. However, the early success of model atmosphere 
codes transformed them into standard tools, and only in the past 
decade have these codes moved beyond simple plane-parallel, 
local-thermodynamic-equilibrium (LTE) models to full three- 
dimensional, statistical-equilibrium codes that can model non- 
LTE physics as well as stellar convection. Unfortunately, com- 
puting power is still limited for calculating large-scale model 
atmosphere grids varying stellar gravity, effective temperature, 
stellar mass and composition. 

A step toward more realistic geometry is achieved by shifting 
from one-dimensional plane-parallel model stellar atmosphere 
codes to one-dimensional spherically symmetric codes, which 
can be used to compute large grids of models atmospheres that 
include physics that is more appropriate to stars where the depth 
of the stellar photosphere is a significant fraction of the stel- 
lar radius, such as evolved giant and supergiant stars and pre- 
main sequence stars. One such code for modeling atmospheres 
assuming spherically sy mmetric geometry is the SAtlas code 
dLester & N eilson 2008). This code is based on the Atlas code 
developed by iKuruca (Q,979), and continues its assumption of 
local thermodynamic and hydrostatic equilibrium. However, the 
radiativ e transfer is com puted assuming spherical geomet ry us- 
ing the iRvbickil dl97lh version of the iFeautrierl (fl964) ray- 
tracing method, while radiative and co nvecting equilibrium is 
enforced using an updated version of the lAvrett & Krookl (1 1963b 
temperature correction method. Models computed using this 
code have been c ompared to spherically-symmetric Phoenix and 
MARCS models dHauschildt et alJl999tlGustafsson et alJ 2008) 
and shown to produce similar results dLester & Neil son 2008; 
iNeilson & Lesteill2008l) . 

In this work we use the grid o f spherical model atmo- 
spheres from Neils on & Lester! (1201 ll) . extended in mass up to 
M — 20 M e . The grid assumes solar composition and spans the 
gravities from logg = -1 to logg = 3 in steps of 0.25, effec- 
tive temperatures from T e g = 3000 to 8000 K and masses from 
M = 2.5 to 20 M Q in steps of 2.5 M and includes models with 
masses M = 0.5 and 1 M . Surface intensities are computed for 
each model at 1000 equally spaced values of p = cos 6, where 8 
is the angle between the vertical direction and the direction to- 
ward a distant observer. Limb-darkening pr ofiles are computed 
for Johnson-Cousins BV7?/7/A'-wavebands dJohnson & Morgan 
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Fig. 2. The limb-darkening coefficient u, used in Eq.Q] applied to 
the Kepler photometric band. Red crosses are the plane-parallel 
model stellar atmospheres, and the blue squares are the spherical 
models. 



ll953 HBesselll2005l) along wi th the CoRot dAuvergne et al.l2 009) 
and Kepler ( Koch et al. 2004) wavebands. Angular diameter cor- 
rections for interferometric observations, gravity-darkening co- 
efficients and various limb-darkening relations are computed us- 
ing these wavelength-integrated intensity profiles. 

3. Limb-Darkening Laws 

An understanding of stellar limb darkening is required to 
model the properties of interferometric, eclipsing binary -star, 
microlensing, and planetary-transit observations. As these ob- 
servations become more precise and more accurate, models of 
stellar limb darkening must also improve. Limb darkening is 
typical ly treated as a simple parametrizatio n as a function of 
(e.g. iFouque et al.ll2010t ICroll et al.ll201 ll) . which makes fit- 
ting the stellar intensity profile much simpler and reduces the 
number of free parameters. The m ost common parametrization s 
are linear and quadratic relations dAl-Naimivlll978i Ivan Ha mme 
H993t lD iaz-Cordoves et al. but other suggested relations 

include a four-parameter relation (IClar et 2000a), a square-root 
relation dWade & Rucinski 1985) as well as exponential and log- 
arithmic relations dClare t 2000a; Claret & Hausc hildt 2003J). 

3.1. Best-Fit Limb-Darkening Laws 

We fit the following limb-darkening relations to the grids of 
plane-parallel and spherically symmetric model stellar atmo- 
spheres: 
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Fig. 1. Kepler-band model intensity profiles (black-solid) predicted for both plane-parallel (left) and spherically symmetric (right) 
model stellar atmospheres with T e ff = 5000 K, \ogg = 2 and M — 10 M . Along with the intensity profiles, best-fit linear (green- 
dashed), quadratic (orange-short-dashed), square-root (blue-dotted), four-parameter (violet-long-dash-dotted), logarithmic (brown- 
short-dash-dotted), and exponential (grey-double-dash) limb-darkening relations are plotted. Bottom panels show the difference, 
A = /model - /law, between model intensities and best-fit limb-darkening laws. 



T/ I(JI \, = 1 - £(1 - J") - h -r~T, Exponential, (5) 

I(fi — 1) 1 - ef 

— — — = 1 - mil — Li) -nu\nu Logarithmic. (6) 

IQx = 1) 

We derive the best-fit coefficients for each of the limb- 
darkening laws using a general least-squares algorithm. This was 
done using the computed surface intensities for the BVRIHK- 
and CoRot- and Keple r-wavebands. Fig. Q] shows the Kepler- 
band intensity profile and corresponding best-fit limb-darkening 
laws for both spherical and plane-parallel model atmospheres 
with the properties T eff = 5000 K, logg = 2 and M = 10 M 
(mass is defined for the spherical model only). The chosen limb- 
darkening laws all fit the plane-parallel model intensity profiles 
well. This is not surprising because plane-parallel model atmo- 
sphere intensity profiles do not deviate significantly from being 
linear, and a linear term is included in all of the chosen limb- 
darkening laws. However, spherically symmetric model stellar 
atmospheres have intensity profiles that are significantly non- 
linear, and the best-fit limb-darkening relations for these inten- 
sity profiles match less well than for the plane-parallel models 
because of this non-linearity. For the model shown in Fig. [T] 
limb-darkening laws predict intensities that vary by A = / mo dei _ 
/law = 0.15 for the spherical model while A < 0.04 for the plane- 
parallel model. Although limb-darkening laws fit the intensities 
of plane-parallel model atmospheres better than spherically sym- 
metric models, the spherical models are more physically realis- 
tic, making them the more appropriate choice to use in modeling 
observations. We explore the uncertainty of the limb-darkening 
fits later. 

We present in Figs. [2] [3] |4] and [5] the coefficients derived 
by least-squares fitting for the limb-darkening laws given by 
Eqs. 1-6 respectively for the Kepler photometric band as a 



function of effective temperature for both plane-parallel and 
spherically symmetric model stellar atmospheres. It is clear that 
more realistic spherically symmetric model stellar atmospheres 
predict limb-darkening coefficients that vary much more as a 
function of effective temperature than those for plane-parallel 
models. For the simplest case of the linear limb-darkening law, 
the M-coefficient determined from plane-parallel models in the 
Kepler-band vary from u - 0.2 to 0.5, whereas spherical mod- 
els with the same effective temperatures and gravities vary 
from u = 0.6 to 1.4. The coefficients predicted for almost all 
limb-darkening laws examined here show the same behavior as 
the limb-darkening coe fficients predicted for a flux-c onserving 
linear+square-root law (Neils on & Lesterl201 lll20L2l) . This uni- 
form dependence of the coefficients on is surprising and sug- 
gests all of these laws carry essentially the same information 
regarding the moments of the intensity and the atmospheric ex- 
tension ab out the stellar atmosphere in question. The one excep- 
tion is the IClarea (120 00a) four-parameterlimb-darkening law, for 
which the coefficients appear to vary much more as a function of 
effective temperature. 

To explore the interdependence of the coefficients, we plot 
in Fig.|6]the Kepler-band ^-coefficient from the quadratic law as 
a function of the a-coefficient. This plot is typical of all the two- 
parameter limb-darkening laws consid ered in this work as well 
as the limb-darkening law employed bv iNeilson & Lesterl(l201 ll 
2012), including the apparent hook in the correlation between 
coefficients. Fig. [6] also plots the values of fa + fa as a function 
of fa + fa for the four-parameter law, again for the Kepler pho- 
tometric band. The correlation for both plane-parallel and spher- 
ical models is readily apparent. A best-fit linear relation to the 
coefficients for spherical models is 

/2,Kepler + AlCepler = -0.989(/i,Kepler + AlCepler) + 1 -05 1 . (7) 
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Fig. 3. Limb-darkening coefficients a and b used in Eq. [2] (left panel), and the coefficients c and d used in Eq. [3] (right panel), all 
applied to the Kepler photometric band. The symbols have the same meanings as in Pig. El 




Fig. 4. Limb-darkening coefficients g and h used in Eq. [5] (left panel), and the coefficients m and n used in Eq. [6] (right panel), all 
applied to the Kepler photometric band. The symbols have the same meanings as in Fig. [2] 



The correlation is different for plane-parallel models for which 
the slope is -0.978 and the intercept is 0.493. 

These correlations are caused by the limb-darkening coeffi- 
cients being linear combinations of various angular moments of 
the intensity. For instance, in plane-parallel model atmospheres 
the moments J = J I(n)dji and K = J I(/j)/x 2 dp are related such 



that J — 3K (Mihalas 1978). In spherical symmetry, this ratio 
varies, causing the moments of the intensity to differ in spheri- 
cal symmetry from those predicted moments for plane-parallel 
model stellar atmospheres. This difference in geometry is re- 
flected in the difference between the zero-points of the relation 
Eq. [7] for spherical models and that for plane-parallel models. 
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Fig. 5. Limb-darkening coefficients f\, fy, fa and f\ used in the IClarel (l2000al) four-parameter law, Eq. 2] applied to the Kepler 
photometric band. The symbols are the same as in Fig. [2] 
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Fig. 6. Correlation between the Kepler-band limb-darkening coefficients for the quadratic law (left) and for the four-parameter law, 
H + /4 as a function of f\ + (right). Red crosses represent coefficients computed from plane-parallel models and blue squares 
spherical models. 



One can potentially use this difference to test observations and cooler stars. We compute the relative error of the limb-darkening 
test model geometry. fit, A, using the relation 



3.2. Error Analysis 

Various limb-darkening laws, such as those given in Eqs. [TJ- 
[6] are fit to the surface intensities computed with model 
stellar atmospheres, and it is important to understand how 
well these laws represent the actual intensities. For instance, 
iDiaz-Cordoves et all (1 19951) argued that a square-root law fit in- 
tensity profiles for hotter stars (T e ff > 9000 K) better than a 
quadratic law, whereas no limb-darkening law is preferred for 



Hf/modelO") ~ht(fj)] 2 



(8) 



which quantifies the deviation of the best-fit limb-darkening law 
from the surface intensities of the model atmosphere. We com- 
pute the relative error for each bandpass as a function of the fun- 
damental stellar parameters for both plane-parallel and spherical 
geometries, and show in Fig. [7] the relative errors as a function 
of effective temperature for fits in the Kepler-band. The relative 
error of the fits for spherical models is greater than the error 
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for plane-parallel fits for all the limb-darkening laws. The errors 
are similar only for T e ff ~ 3500 K, where the spherical model 
atmospheres predict intensity profiles that are closest to being 
linear, with the error of the linear limb-darkening law approach- 
ing a minimum value. This result appears to suggest that these 
limb-darkening laws are inappropriate for fitting light curves and 
interferometric observations, but this is not true. 

There are a number of issues with how the relative error is 
computed and what the error tells us, such as how the limb- 
darkening laws are defined, how they are fit to the surface in- 
tensities and the effect of sampling. 

- Defining limb-darkening laws: The intensity profiles com- 
puted using the plane-parallel and spherical model atmo- 
spheres employed in this work are normalized with respect 
to the central intensity so that = 1) = 1. Furthermore, all 
limb-darkening laws, except the exponential law, are defined 
so that the I{fi = 1) = 1, regardless of the values of the best- 
fit coefficients. As a result, every fit to an intensity profile is 
anchored to the center of the stellar disk before representing 
the remainder of the intensity profile. This definition alone 
results in a perfect fit to the center of the stellar disk and a 
deteriorating fit as — > as the intensity profile deviates 
from the assumed structure of a particular limb-darkening 
law. 

- Fitting limb-darkening laws: Limb-darkening laws are typ- 
ically fit to intensity pro files using a least-square algorithm. 
iNeilson & Lesterl d201 ll) showed that the best-fit coefficients 
for a given law are functions of weighted integrals of the 
intensity profile. For example, the linear limb-darkening co- 
efficient from Eq. 1 is a function of the mean intensity, J, and 
the stellar flux, "H = J I(jx)fidfx, and both of these quantities 
are more sensitive to the central intensity than to the much 
smaller intensity near the limb. As with the definition of the 
limb-darkening laws, using a least square fitting algorithm 
fits the central part of the intensity structure better. Similarly, 
one might fit limb-darkening coefficients by enforcing flux 
conservation, but because the flux is the //-weighted integral 
of the intensity, any flux-conserving fit is constrained weakly 
by the intensity at the stellar limb relative to the intensity 
near the center of the stellar disk. 

- Sampling Issues: Sampling is the most important of the 
three issues affecting the co mputed error of th e fit of the 
intensity profile. F or instance, IWade & Rucinskil (1 19851) and 
Heyrovsky (2007) noted that fitting an intensity profile that 
is uniformly sampled in // has a larger error than fitting 
the same profile that is uniformly sampled in r — sin 6 — 
sin(cos _1 //). Uniform r-spacing emphasizes the intensity 
profile near the center of the disk while a uniform //-spacing 
emphasizes the limb. Adopting any of the limb-darkening 
laws presented here, that law will fit the central part of the 
stellar surface more precisely than the limb because of the 
normalization at the center of the disk. If, in addition, the 
surface intensity is sampled uniformly in r, that will give 
added weight to the central region. These two factor com- 
bine to make the computed error of the fit smaller. Similarly, 
iHowarthl (1201 lbl) found that limb-darkening coefficients de- 
rived from planetary transits with large impact factors do 
not agree with model stellar atmosphere predictions. This is 
because the planet passes across only the limb of the star 
and not the c enter, therefore prob ing only part of the inten- 
sity profiles. IClareti (2008, 2009) also found disagreement 
between theoretical limb-darkening coefficients and empiri- 
cal coefficients measured from eclipsing binary light curves 



and comparisons to the planetary system HD 209458. Limb- 
darkening coefficients from stellar atmosphere models fit the 
whole profile yielding different results. 

The combination of these three factors lead to calculated er- 
rors that are relative and not an absolute measure of the quality of 
the fit. In this work, differences in the error between fits to plane- 
parallel and spherically symmetric model stellar atmosphere in- 
tensity profiles computed with the same properties are due solely 
to differences in the intensity profile near the limb where the 
spherical models provide more realistic predictions. Therefore, 
the error analysis suggests that the various limb-darkening laws 
lack the necessary complexity to precisely fit in tensity profiles 
from spherical models. The only exception is the IClareti 



pro tiles 

(2000a) 



four-parameter law, which fits the laws best, but appears to have 
unique properties. 

4. Gravity Darkening Coefficients 

Rapid rotation distorts the shape of a star, making it aspheric, 
with flattened pol es and a bulged equator. As shown first by 
Ivon ZeipeH d 19241) . the gravity and effective temperature vary 
in a coordinated way across the stellar surface such that at 
any point the effective temperature is proportional to the ef- 
fective gravity, T e g ~ gf^ 4 , where fi\ — 1 for radiative stars. 
However, this value of fi\ is valid only for bolometric radia- 
tion, and iKopall d!959h lat er derived monochromatic gravity- 
darken ing corrections, y(X). Claret J 2000al) . lciaret & H auschildt 
(2003}) and lClaret & Bloemenld201 ll) have computed waveband- 
dependent gravity-darkening corrections as a a function of the 
central intensity of the star, as well as the gravity, effective 
temperature and the variable, fi\ from plane-parallel models. 
iBloemen etail d2011l) derived 



y(X) 



d In 1(A) 
Sing 



d\nT e s\ld\nI(A) 



ding j\d\nT e ff 



(9) 



and noted that (d In T e s/d Ing) = fl\ /4. The variable fl\ is a func- 
tion of effective temperature, but for the purpose of this analysis 
we assume ji\ = 0.2 for r e ff < 7500 K and fi\ — 1 for hotter 
stars. However, the value of fi\ based on von Zeipel's theorem 
is not strictly valid for radiative or convective stellar envelo pes 
dClaretll2000bt lEspinosa Lara & Rieutordll201 U IClaretll2012l) . 

In Fig. [8] we plot the V-band values of each intensity deriva- 
tive for each model stellar atmosphere in Eq. [9] as well as y(A) 
computed for the assumed values of f}\. We find that plane- 
parallel and spherically symmetric model stellar atmospheres 
predict similar gravity-darkening coefficients for T^g > 4000 K, 
but there are significant differences for cooler stars. We interpret 
these differences for the cooler stars as consequences of both 
surface convection and the shift from the negative hydrogen ion 
to titanium oxide as the dominant opacity source. Both plane- 
parallel and spherical model intensities show greater variation 
at these cool effective temperatures, but the intensity profiles of 
spherically symmetric model atmospheres vary more than that 
of plane-parallel model atmospheres. 

While the most significant differences between spherical and 
planar model predictions of gravity-darkening coefficients are 
at lower temperatures, the gravity-darkening coefficients com- 
puted from spherically symmetric models are greater than those 
of plane-parallel models for every effective temperature. For ex- 
ample, a spherically symmetric model with T e g = 8000 K has a 
V-band gravity-darkening coefficient of yy - 0.165 while the 
plane-parallel model with the same effective temperature has 
yv — 0.14. The difference is small but systematic. 
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Fig. 7. The error of the best-fit limb-darkening relation, defined by Eq. [8] for every model atmosphere (red crosses represent plane- 
parallel models, blue squares spherical models) for each of the six limb-darkening laws at Kepler-band wavelengths. 



5. Angular Diameter Corrections 

Interferometric observations measure the angular diameter of 
a star along with its limb-darkening profile, but, unfortunately, 
the measured angular diameter and limb-darkening profiles are 
not independent quantities. This is especially tru e when the 
measu red visibilities do not probe the second lobe. iDavis et ail 
(2000) measured stellar angular diameters from interferomet- 
ric observations by assuming that the stellar intensity profile 
is uniform, i.e. the intensity at any point on a stellar disk is 
equal to the central intensity. In that case, the uniform-disk an- 
gular diameter can be directly fit to the observed visibilities 
and then converted to a lim b-darkened a ngular diameter using 
model stellar atmospheres. IDavis et all (l2000|) computed cor- 
rections using plane-parallel Atlas models ( Kuruczlll993l) and 
found k = 6>ud/#ld = 0.91 to 0.98 in the wavelength range 
A = 400-800 nm. These limb-darken ing corrections have been 
applie d to observations of Cepheids (Galle nne et al.ll2012h and 
Sirius dDavis et al.ll201 ll) for example. 



We comp ute angular diameter corrections using the recipe 
described by Marengo et al. (2004), where we assume a limb- 
darkened angular diameter of #ld = 1 mas to compute interfer- 
ometric visibilities from a model atmosphere intensity profile. 
That synthetic visibility is then fit by a uniform-disk angular di- 
ameter. The best-fit uniform-disk angular diameter is then equiv- 
alent to the theoretical angular diameter correction. We compute 
angular diameter corrections for the Johnson-Cousins BVRIHK 
wavebands and show the corrections for the V- and A'-bands in 
Fig. [9] as a function of effective temperature for plane-parallel 
and spherically symmetric models. Corrections from spherical 
models clearly differ from corrections from plane-parallel model 
atmospheres. Intensity profiles from plane-parallel model stel- 
lar atmospheres predict corrections in the narrow range from 
k = 0.97 - 0.99 in V-band and approaches unity for longer wave- 
lengths. We show in Fig. [9] the V- and /f-band angular diame- 
ter corrections as function of effective temperature and gravity 
for plane-parallel and spherically symmetric model stellar atmo- 
spheres. 
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Fig. 8. V-band central intensity derivatives and gravity-darkening coefficients as function of effective temperature (left) and gravity 
(right) computed from plane-parallel (red crosses) and spherically symmetric (blue squares) model stellar atmospheres. 



Fits to spherically symmetric model atmospheres suggest 
significantly different angular diameter corrections as functions 
of both effective temperature and gravity. The V-band correc- 
tions from spherical models, denoted k s , range from k s = 0.67 
to 0.95, with no overlap with the plane-parallel model predic- 
tions. The /T-band corrections show similar behaviors except 
that spherical and planar corrections overlap somewhat. These 
results suggest that using plane-parallel model atmosphere cor- 
rections systematic ally underestimates the st ellar angular diame- 
ter. For instance, Mozur kewich et all d2003) presented uniform- 
disk angular diameters for a sample of 85 stars, along with limb- 



darkened angula r diameters correc ted u sing limb-darkening 
coeffic ients from IClaret etall (1 1995b and iDiaz-Cordoves et all 
(Il995h . Their angular diameter corrections vary from k = 0.89 
to a: 1, consistent with the values found here for plane-parallel 
model atmospheres. 

Of particular interest are the results of Moz urkewich et all 
d200l for a Persei (F5Ib), for which they measured T e s = 
6750 K, and for e Ge minorum (G8 lb), w h ich w as measured to 
have r e ff = 4485 K. iMozurkewich et aD (|2003) measured the 
uniform-disk angular diameters at 550 nm to be 2.986 + 0.042 
for a Per and 4.467 + 0.115 mas for e Gem. Using these they 



8 



H. R. Neilson: Limb darkening in red giant stars 




CD 

Z> 
CD 

T3 
£= 
CO 

CD 



CD 

Z> 
CD 

T3 
c= 

CO 
J2 



1 

0.95 
0.9 

0.85 
0.8 

0.75 

0.7 

0.65 
1 

0.95 
0.9 

0.85 
0.8 

0.75 
0.7 

0.65 



i 1 1 1 1 r 



^^j*tt 




3000 4000 5000 6000 7000 8000 
T ef f (K) 



-2-10 1 2 3 
Log g (cgs) 



Fig. 9. Interferometric angular diameter correction computed in V-band (top) and Zf-band as functions of effective temperature 
(left) and gravity (right). Corrections computed from plane-parallel model atmospheres are denoted with red x's and spherically 
symmetric models blue squares. 



computed limb-darkened angular diameters of 3 . 1 88 ± 0.035 and 
4.703 ±0.047 mas, respectively. Our spherically-symmetric mod- 
els with logg = 1.5 and M = 10 M yield V-band angular- 
diameter corrections of 0.929 for a Per and 0.916 for e Gem. 
Applying these to the uniform disk measurements gives larger 
limb-darkened angular diameters: $ld = 3.214 mas for a Per 
and 6>ld = 4.877 mas for e Gem. The spherical correction 
for a Per yields a value for #ld that is marginally consistent 
with the angular diameter found using plane-parallel correction, 
wh ereas the limb-darkened an gular diameter of e Gem measured 
bv iMozurkewich et al.1 (2003) is almost 4% smaller than what 
would be predicted by applying spherical model corrections. 
This difference may appear to be small but this underestimate 
is systematic. 

As a test, we check how the angular diameter corrections 
vary as function of stellar mass. Because models with low effec- 
tive temperature but relatively high gravity appear to predict the 
smallest corrections, we hold T e s = 3500 K and logg = 2. The 
angular diameter corrections are shown in Fig. [I0]as a function 
of stellar mass for the six Johnson-Cousins wavebands consid- 
ered in this work. The figure suggests that the corrections are 
insensitive to the mass of the stellar model except for low-mass 
(M < 1 M Q ) models. This is reassuring and suggests that when 
applying these corrections, one can ignore the stellar mass. The 



difference between limb-darkening profiles and angular diame- 
ter corrections is s mall and consistent with previous results by 
iLester et al.l(l2013l) . 

6. Summary 

In this work, we present model atmosphere intensity profiles for 
the BVRIHK, CoRot and Kepler passbands from both plane- 
parallel and sph erically symmetric geometries based on mod- 
els computed bv lNeilson & Lester! d2011l |2012|) . We fit a num- 
ber of limb-darkening laws to these intensity profiles, as well 
as compute gravity-darkening coefficients and angular diameter 
corrections for interferometry. We test how these fits vary as a 
function of model atmosphere geometry and compile tables of 
limb-darkening coefficients, gravity-darkening coefficients and 
angular diameter corrections that can be applied to observations. 

We consider six limb-darkening laws in this work: linear, 
quadratic, square-root, four-parameter, exponential and logarith- 
mic. These laws fit the intensity profiles from plane-parallel 
model atmospheres well, but not the intensity profiles of the 
spherical mode ls bas e d on co mputed relative errors. The one ex- 
ception is the IClared (12000a) four-parameter law, for which the 
difference between the spherical model intensity profiles and the 
predictions of the fitting law is small enough to still be applicable 
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Fig. 10. Interferometric angular diameter corrections as a func- 
tion of waveband and stellar mass for spherically symmetric 
model stellar atmospheres with log g - 2 and T e g = 3500 K. 



Table 1. Summary of limb-darkening coefficient, gravity- 
darkening coefficient and interferometric angular diameter cor- 
rection tables found online. 



Name Geometry 



Typf 

Table2 Spherical Linear Limb Darkening Eq. [1J 

Table3 Spherical Quadratic Limb Darkening Eq. [2] 

Table4 Spherical Square Root Limb Darkening Eq. [3] 

Table5 Spherical Four-parameter Limb Darkening Eq.|4] 

Table6 Spherical Exponential Limb Darkening Eq. f5] 

Table7 Spherical Logarithmic Limb Darkening Eq. [6] 

Table8 Planar Linear Limb Darkening Eq. [TJ 

Table9 Planar Quadratic Limb Darkening Eq. [2] 

Table 10 Planar Square Root Limb Darkening Eq. [3] 

Table 11 Planar Four-parameter Limb Darkening Eq.|4] 

Table 12 Planar Exponential Limb Darkening Eq. f5] 

Table 13 Planar Logarithmic Limb Darkening Eq. [6] 

Table 14 Spherical Gravity Darkening 

Table 15 Planar Gravity Darkening 

Tablel6 Spherical Angular Diameter Corrections 

Table 17 Planar Angular Diameter Corrections 

Note 1. Tables listed here can be retrieved electronically from the CDS. 



to observations, although the law still fits the spherical profiles 
more poorly than the plane-parallel intensities. 

While those predicted errors are useful for comparing fits to 
planar and spherical model intensity profiles, they are not ideal 
for studies of actual limb darkening. Best-fit limb-darkening co- 
efficients depend on the definition of the laws, all of which an- 
chor the fit to = 1) = 1, making the fit sensitive to the 
sampling of the intensity profile as well as to the method for 
fitting the data. Because intensity profiles for spherical models 
are more complex, the fitting error is greater than the error for 
simpler plane-parallel model intensity profiles. However, spher- 
ically symmetric model atmospheres are a more realistic repre- 
sentation of actual stellar atmospheres, meaning they are better 
suited for limb darkening studies. 

Fits to the four-parameter limb-darkening law also show cor- 
relations between the limb-darkening coefficients; we find that 
the linear combination of the four coefficients are approximately 
constant, with that constant being a function of the atmosphere's 
geometry. This result suggests that the linear combination of the 
observed coefficients for the four-parameter law provides a sim- 
ple test of whether the observations are probing the edge of the 
stellar disk, i.e. sphericity. 

We also predict wav elength-dependent gravity- darkening co- 
efficients based of the IClaret & BloemerJ (1201 ll) prescription. 
Unlike the limb-darkening coefficients, the gravity-darkening 
coefficients are less dependent on model atmosphere geometry. 
This is because the gravity-darkening coefficients depend on the 
change of the central intensity with respect to effective temper- 
ature and gravity, hence the difference between atmospheres for 
the same geometry. Gravity-darkening is also a function of the 
central intensity, which is insensitive to model geometry. The 
spherically symmetric gravity-darkening coefficients are similar 
to plane-parallel coefficients for T e g > 5000 K and begin to di- 
verge for cooler stellar atmosphere models. Only at the coolest 
effective temperatures, 3000 K< T e g < 4000 K, is the geometry 
of the model atmosphere important, with the spherically sym- 
metric coefficients being approximately an order-of-magnitude 
greater than those predicted from plane-parallel model atmo- 
spheres. 

Unlike the gravity darkening coefficients, the interferomet- 
ric angular-diameter corrections do depend on geometry. For 



plane-parallel model atmospheres the angular-diameter correc- 
tions vary from about 0.95 - 1, whereas the corrections for 
spherically symmetric model atmospheres vary from 0.67 - 1. 
Previous analyses had assumed that corrections from plane- 
parallel models are applicable to all stars, but this is not true. At 
low gravity, log g < 3, spherically symmetric corrections deviate 
significantly from plane-parallel model predictions. The differ- 
ence between spherical and plane-parallel models is a function 
of both gravity and effective temperature and also appears to vary 
as a function of stellar mass. 

The angular-diameter corrections, limb-darkening and 
gravity-darkening coefficients are publicly available as online ta- 
bles. Each table has the format T ef! (K) , logg, M (M ), R (R Q ) 
and L (L Q ) and then the appropriate variables for each waveband, 
such as linear limb-darkening coefficients. Tables of gravity- 
darkening coefficients also contain values of the intensity deriva- 
tives with respect to gravity and effective temperature. For plane- 
parallel models, values of mass, radius and luminosity are pre- 
sented as zero in the tables. We list the properties of these tables 
in Tab.Q] that are archived in electronic form at the CDS. Model 
atmosphere intensity profiles are also archived at the CDS 

Techniques such as optical interferometry, microlensing ob- 
servations, planetary transit and eclipsing binary observations 
are continuously improving the measurements of stellar limb 
darkening needed to test model stellar atmospheres and the 
physics assumed in their calculation. At lower gravities, these 
observations require the more physically realistic spherically 
symmetric models to constrain stellar properties. The predicted 
limb-darkening coefficients, gravity-darkening coefficients and 
angular diameter corrections from spherically symmetric S Atlas 
models are new tools that for aiding analyses of these observa- 
tions. 
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